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Abstract 

In previous work [L. Blanchet and A. Le Tiec, Phys. Rev. D 78, 024031 (2008)], a model of 
dark matter and dark energy based on the concept of gravitational polarization was investigated. 
Lhis model was shown to recover the concordance cosmological scenario (A-CDM) at cosmological 
scales, and the phenomenology of the modified Newtonian dynamics (MOND) at galactic scales. 
In this article we prove that the model can be formulated with a simple and physically meaningful 
matter action in general relativity. We also provide alternative derivations of the main results of 
the model, and some details on the variation of the action. 
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I. INTRODUCTION AND MOTIVATION 



The concordance cosmological model A-CDM brilliantly accounts for the precise measure- 
ments of the anisotropies in the cosmic microwave background (CMB) pQ, for the formation 
and growth of large scale structure as seen in deep redshift [2] and weak lensing [3] surveys, 
and for the fainting of the light curves of very distant supernovae jU [5]. The paramount 
conclusion is that the total mass-energy content of the Universe is made by f^b — 4% of ordi- 
nary (essentially baryonic) matter, fl^m — 23% of cold dark matter (CDM), and fide — 73% 
of dark energy which could be in the form of a cosmological constant A. However, no exper- 
imental claim of direct discovery of a CDM particle has been confirmed, and the attempts 
at interpreting A in terms of fundamental quantum mechanics have failed. 

With the advent of high precision cosmic V-body simulations (see [BJ for a review), the 
cosmological model has been extrapolated to the smaller scale of galactic systems, and sug- 
gests the existence of a specific CDM density profile around galaxies [7]. However, the 
simulated CDM halos face severe challenges when compared to observations. Most prob- 
lematic is the generic formation of cusps of dark matter in the central regions of galaxies, 
while the rotation curves seem to favor a constant density profile in the core [H [9]. We 
mention also the prediction of numerous but unseen satellites of large galaxies [101 [11], and 
the recent evidence [12] that tidal dwarf galaxies are dominated by dark matter — a fact 
which is at odds with the CDM tenets [13]. Furthermore, the most important challenge is 
that the CDM scenario falls short in explaining in a natural way Milgrom's law [T4"l [THJ HE] , 
namely that the need for dark matter arises only in regions where the typical acceleration 
of ordinary matter (or, equivalently, the typical value of the gravitational field) is below 
some universal constant acceleration scale a® ~ 1.2 x 1CT 10 m/s 2 . This law manifests itself 
particularly in the flat rotation curves of spiral galaxies, and in the baryonic Tully-Fisher 
relation. No convincing mechanism for incorporating an acceleration scale such as ao in the 
V-body simulated CDM halos has been found. Although it is possible that some of these 
problems will be solved within the CDM approach [TTJ [18], it is very important to consider 
alternative solutions. 

The most successful alternative approach to the problem of dark matter in galactic halos 
is MOND — Milgrom's modified Newtonian dynamics [HI [TJ)J [IB], which insists that there 
is no dark matter and we instead witness a violation of the Newtonian law of gravity. In 
MOND the true gravitational field g experienced by ordinary matter (stars and gas) differs 
from the Newtonian one, and obeys the modified Poisson equation [19] 



V- (fig) = -47rp b . (1.1) 

We use bold-face notation to represent ordinary three-dimensional vectors and pose G — 1. 
Here pb is the density of baryonic matter, and ft is the MOND function which depends 
on the norm g — \g\ of the gravitational field. In the regime of weak gravitational fields, 
g <C a , we have fi(g) = g/a + 0(g 2 ), while fi(g) — > 1 when g ^> ao, so as to recover the 
usual Poisson equation. Various forms of the interpolating function « have been proposed 
to fit observations in the best way 



The ability of the formula (1.1) to reproduce a wide variety of phenomena associated 



with dark matter halos is tremendous (see e.g. [221 123])- However, because (1.1) is non- 
relativistic, it does not allow one to answer questions related to cosmology. In particular, it is 
a great challenge to find a theory reproducing both MOND at galactic scales and A-CDM at 
cosmological scales. A number of relativistic field theories have been proposed, recovering 
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(1.1) in the non-relativistic limit, and sharing with MOND the idea that dark matter is 
an apparent reflection of a fundamental modification of gravity. The prime example of 
such modified gravity theories is the tensor-vector-scalar theory (TeVeS) of Bekenstein and 
Sanders [2U EE]. Interesting connections between TeVeS and the class of Einstein-aether 
theories [27] have been found [2SJ [2JJJ ED] • Modified gravity theories are rather complicated 
extensions of general relativity (GR), and are for the moment not connected to fundamental 
physics. Moreover, they do not account for all the mass discrepancy at the intermediate scale 
of galaxy clusters |31j . To resolve this difficulty a component of hot dark matter (HDM) in 
the form of massive neutrinos has been invoked (32j [33] . At cosmological scales the modified 
gravity theories also have some problems at reproducing the observed CMB spectrum [34], 
even when using a component of HDM. 

The approach we propose below is able to successfully address both cosmological and 
galactic scales. We advocate that a non-standard form of dark matter may exist, while 
keeping the standard law of gravity (GR) unchanged. The physical belief of this alterna- 
tive approach is the striking analogy between MOND and the electrostatics of (non-linear 
isotropic) dielectric media [35] . Indeed, the MOND equation (1.1) can be interpreted as the 
standard Poisson equation if the gravitational field is sourced by baryonic matter and by 
a "digravitational" medium playing the role of dark matter. The density of "polarization 
masses" in this medium is then p po i = — V • (anticipating the notation adopted below), 
where II ^ denotes the polarization field, which must be aligned with the local gravitational 
field, 

n ± = -X&g. (1.2) 

Here x = A* ~ 1 denotes the "gravitational susceptibility" coefficient of the medium, while \i 
can be viewed as a "digravitational" constant. It was argued [35] that in the gravitational 
case the sign of x should be negative, in agreement with what MOND predicts; indeed, we 
have n < 1 in a straightforward interpolation between the MOND and Newtonian regimes, 
hence x < 0. Furthermore, arguments were given showing that the stability of the dipolar 
medium requires the existence of some environment-dependent internal non-gravitational 
force. More precisely, the force has to depend on the polarization field, i.e. the density of 
dipole moments. 

Motivated by the previous interpretation of MOND we present in Section [IT] an action 
principle for dark matter viewed as the gravitational analogue of a polarizable dielectric 
medium. In Section III we show that this model is currently viable since it is in agreement 
with the standard cosmological scenario at large scales and recovers MOND at galactic 
scales. Some details regarding the variation of the action are relagated to Appendix [A] 



II. MODEL OF DIPOLAR DARK MATTER AND DARK ENERGY 

In previous work [35] (hereafter paper I; see also [37J for an earlier attempt) we proposed 
a relativistic model of dark matter and dark energy based on a particular concept of gravi- 
tational polarization. In contrast to modified gravity theories, the model should be viewed 
as a modified matter theory. The idea that the phenomenology of MOND could arise from 
the CDM paradigm has been previously discussed [381 ES] . However here we shall consider 
a true modification of the physics of dark matter, drastically different from CDM (see also 
[40] for an alternative approach in a related spirit). 



3 




FIG. 1: The dipolar fluid is entirely characterized by (i) its four- velocity u M tangent to the world-line x^(t), 
(ii) its rest mass density a, and (iii) the dipole moment The polarization II j_ = cr£j_ is build from the 
norm £j_ of the projection £^ of the dipole moment £ M onto the space-like hypersurface E t orthogonal to . 



In paper I we showed that this particular model of modified dark matter permits recov- 
ering the phenomenology of MOND in a natural way, while being in agreement with the 
cosmological A-CDM model. The aim of this article is to prove that the model can be 
reformulated from a simple and physically meaningful matter action. 

The dipolar medium is described as a fluid with mass current J M = au^, and endowed 
with a dipole moment vector Here = dx^/dr is the time-like four- velocity of the fluid, 
with dr = yj— g^dx^dx" being the proper time (we pose c = 1). The rest mass density 
reads o = ^J—J^Jv, and the mass current is conserved, i.e. 

V M J" = 0, (2.1) 

where denotes the covariant derivative associated with the metric g^. The dipole mo- 
ment £ M has the dimension of a length, so that it is more like a displacement vector; the 
associated polarization field then reads IP = cr£ M . We have in mind that ^ and IP are 
effective variables resulting from an average performed at some macroscopic scale. 

The action of the dipolar dark matter is of the type S = f d 4 x y/—g L. It will be added to 
the Einstein-Hilbert action for gravity, and to the standard actions of all the other matter 
fields (baryons, photons, neutrinos, etc). We find that the Lagrangian consists of three 
terms: a mass term a in the ordinary sense (as for CDM), an interaction term built from 
the coupling between the current J M and the dipolar field and a potential scalar function 
W characterizing an internal force acting on the dipolar particles, and depending on the 
polarization field. It explicitely reads 

L = -o + J^t - W(II ± ) . (2.2) 

Both J M and £ M will be considered as dynamical variables, to be varied independently when 
applying the principle of stationary action (cf. Appendix LAj for details). Here £ M = D£ M /dr = 
u v V v ^ denotes the proper time derivative of the dipole moment. Note that £^ is coupled 
to the current J M like an external field would be. However, here the dipole moment 
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is an internal field. The potential W is assumed to depend solely on the norm of the 
projection of the polarization field perpendicular to the four- velocity. Thus, IT^ = cr£ iwhere 
Ci = a/-L^ C^ 7 ) with the usual orthogonal projector _L^= 57^ + m^m^ (cf. Fig. lit) . As 
discussed in Sec. |TJ the introduction of such an environment-dependent potential is motivated 
by the previous interpretation of the MOND phenomenology resulting from the mechanism 
of gravitational polarization. 



A crucial point is that £ M can be replaced in the coupling term of (2.2) by its orthogonal 
projection £^ = £ u without changing the dynamics. Indeed, a short calculation reveals 
that J M £ M = + V f j,(Jv , u u £ u ), so that, because a pure divergence can be dropped from 
the Lagrangian, the only physical degrees of freedom are the three independent components 
of the vector £j\ which is space-like (we denote £^ = D^/dr). This is to be contrasted with 
TeVeS and Einstein-aether theories which are based on a fundamental time-like vector field. 

To obtain the equation of motion of the dipolar fluid we vary the action with respect to 
the dipole moment variable and get 



u 



= W , (2.3) 



where ii p = Dw^/dr is the four-acceleration, £^ =€±/€± is t ne um t direction along and 
W' = dW/dn^. The motion is non-geodesic because of the internal force density JF M caused 
by the dipole moment £ M . 

The variation with respect to J M yields the equation of evolution for the dipole moment. 



The constraint that the matter current is conserved, Eq. (2.1), is to be satisfied during the 
variation and we apply a convective variational procedure (see Appendix [A| . Defining for 
convenience the "linear momentum" Q p = ^ + (1 + 2£_|_W / ), we obtain 

QH = I y/x (yy _ n ± W') - ftfl^uV • (2.4) 

This tells how the variation of the dipole moment should differ from parallel transport along 
the fluid's worldline. The first term on the right-hand-side (RHS) looks like a pressure 
term, while the second term represents the analogue of the standard coupling to Riemann 
curvature for spinning particles in GR [4Tj |4"2"] . Finally, varying with respect to the metric, 
we get the stress-energy tensor 

T v» = q(hj») - V p ( H p ± u^ - u p Ii [ l u v) ^j - (W - IL ± W) ■ (2.5) 

The RHS is made of a monopolar term associated with Q^, while the second term is (minus) 
the divergence of a "polarization" tensor and is of a dipolar nature. Being proportional 
to the metric, the third term on the RHS will be related to a fluid of dark energy. We 



have Vj/T^" = as a consequence of (2.3)-(2.4). We observe, in agreement with our earlier 



argument at the level of the Lagrangian, that all equations depend in fine only on the 
perpendicular projection £^ = of the dipole moment. 



The equations of motion (2.3) and evolution (2.4), and the stress-energy tensor (2.5), 
turn out to be exactly the same as in the model of paper I [see (2.20)-(2.21) and (2.24) 
there]. Those equations were derived starting from the more complicated Lagrangian given 
by (2.7) in paper I, and sharing some common features with the one for particles with spin 
moving in an arbitrary background [42] . Furthermore, they were obtained after imposing a 



particular choice of solution satisfying some consequence of the initial equations (5 = 1 in 
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the notation of paper I). Despite this rather complicated way to derive them, it was found 
that the equations provide the sensible physics for a successful model of dark matter and 
dark energy. We have now proved that the same equations derive directly (without any 



further assumptions) from the remarkably simple Lagrangian (2.2), which lends itself better 
to physical interpretation. 



III. RECOVERING THE STANDARD COSMOLOGICAL MODEL AND MOND 

We now review the main consequences of this model, presenting alternative versions of 
most arguments compared to paper I. To achieve agreement with MOND and with A-CDM 
(to first-order cosmological perturbations), we have to fine-tune the potential W in the 
action. Indeed, we find that W is "phenomenologically" determined up to third order in 
an expansion when the polarization field 11^ tends to zero. Physically, this corresponds to 
<C ao, which in turn will mean that gravity is weak, g < ao, like in the outskirts of 
a galaxy or in a nearly homogeneous and isotropic cosmology. In this regime W takes the 
anharmonic form 

A 1 fi7T 2 

w(n ± ) = - + 2vr ni + ^- ni + 0(n 4 ± ) . (3.1) 

The minimum is directly related to the cosmological constant A, and the deviations from 
that minimum are fixed by the agreement with MOND; in particular ao parametrizes the 
third-order deviation (see Fig. [2j. 

Let us assume, following paper I, that the theory depends only on one new fundamental 
scale — the constant MOND acceleration a Q . When entering the MOND regime, U±/ao is 
of order one, therefore W naturally scales with a§. If W is to come from some fundamental 



theory, we expect that the dimensionless coefficients in the expansion (3.1) after global 
rescaling by ciq should be of the order of one. In particular, A should itself be of the order 
of Oq. As is well known [33], the current astrophysical measurements verify the "cosmic 
coincidence" that A ~ Oq. This is a natural consequence of our model. 



A. First-order cosmological perturbations 

We now turn to the application at early cosmological time, where we consider a linear 
perturbation around an homogeneous and isotropic Friedman-Lemaitre-Robertson- Walker 
(FLRW) universe. Since the dipole moment £^ is space-like, it will break the spatial isotropy 
of the FLRW background, and must necessarily belong to the first-order perturbation, which 



we indicate by ^ = 0(1). For instance, from (3.1) we find that the internal force is also 
of first order, JP 1 = 47rlF[ + 0(2). At that order the stress-energy tensor (2.5) simplifies 
very much, and can be decomposed into dark energy and dark matter components, namely 
rpnu _ _|_ jim^ w here the dark energy is simply given by the cosmological constant, 
T£ = -^g^ + 0(2), while the dark matter reads 

TZ = pW + 0(2). (3.2) 

Here p = a — V M IT^ is the energy density of the dark matter fluid, and tt M = u^+^ — ^S/yU^ 
(i.e. tt M = m m — w^, where is the Lie derivative) is an effective four-velocity field, 
which satisfies u^u^ = — 1 + 0(2) and the approximate conservation law V At (ptt M ) = 0(2). 
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Newtonian regime 
/ 

/ 




FIG. 2: The potential W is function of the polarization field II j_. It has an anharmonic behavior in the 
MOND regime II_l <C ao, and its minimum is given by the cosmological constant A. The leading order 



non-harmonicity is parametrized by the MOND acceleration scale ao [cf. Eq. (3.1 1]. The Newtonian regime 
II_l ^> ao is discussed in details in paper I. 



This shows that, at linear order, the dark matter cannot be distinguished from a pressureless 
perfect fluid; in particular the fluid's motion is geodesic, u u V v u^ = 0(2). Therefore, the 
model makes the same predictions as the A-CDM cosmological model at linear order (see 
paper I for more details). In particular, adjusting the background value of p (namely p 
such that p = p + 0(1); notice that p — a) to the measured value of dark matter today, 
Q dm ~ 0.23, and choosing A in such a way that the dark energy contribution is fl^e — 0.73, 
we are in agreement with the observed fluctuations of the CMB. To be more precise, the 
linearized perturbation equations, given by (3.48)-(3.49) in paper I, are identical with those 
of A-CDM with no additional degrees of freedom, since the dipole moment has been absorbed 
at linear order into the effective vector field u fl and mass density p. Therefore the model 
reproduces both the location and the height of the peaks of the CMB. 

At non-linear order in cosmological perturbations, the model should differ from the stan- 
dard A-CDM scenario. The fluid's dynamics will no longer be geodesic. We expect that the 
formation of large scale structures will be triggered not by the monopolar part a of dark 
matter, which should not cluster much (see below), but by the dipolar component present 
in p = a — V^n^. This should be checked performing numerical simulations in cosmology. 



B. Non-relativistic limit 



Next we investigate the dipolar medium around a typical galaxy at low redshift. For 
this purpose it is sufficient to consider the non-relativistic (NR) limit of the model, when 
c — ► +oo. Working at the level of the Lagrangian (2.2), we find that £ = sf—g L reduces in 
this limit to 



c 



NR 



w(n_, 



(3.3) 
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Here v is the fluid's ordinary velocity, £j_ is the dipole moment vector, and g = VZ7 is the 
Newtonian gravitational field with U the gravitational potential. The rest mass density cr* 
satisfies the usual continuity equation d t a* + V • J* = 0, where the current reads J* = a*v. 
We denote by d/dt = dt+v-'V the convective time derivative, so that e.g. dcr*/dt = — <t*V-i>. 
We discarded for convenience the rest mass term (— a*) in the NR Lagrangian (3.3 ). From the 



coupling term in the Lagrangian (2.2 ), we recover in the NR limit (3.3 ) the gravitational 
analogue g ■ U± of the coupling of the polarization field to an external electric field. Notice 
also the extra term J+ ■ d$,±/dt in (3.3), which arises in the gravitational case. 

The NR equations of motion and evolution are obtained by varying the Lagrangian (3.3) 
with respect to £ and £° (note that £ 



dv 

dt 
d 2 ^ 
dt 2 
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1 



£ — v£°), and J*. We get 



v(w-n ± W) + (£l-v)0 



(3.4) 
(3.5) 



in agreement with the NR limits of (2.3) and (2.4). The gravitational equation follows from 
adding to (3.3) the Newtonian Lagrangian Cjj = — ^-V£7- VZ7 (coming from the NR limit 



of the Einstein-Hilbert action in GR) and the contribution of baryons. Varying with respect 
to U gives 

V-(0-47rnL) = -47r(pb + <T*). (3.6) 

We proposed in paper I a mechanism by which the dipolar medium does not cluster 
as much as baryonic matter during the cosmological evolution. This is supported by an 
exact solution of (3.4)-(3.6), valid in spherical symmetry, where the dipolar fluid has zero 
velocity, v = 0, and a constant mass density a* (see Appendix A in paper I). The dipole 
moments remain at rest because the gravitational field g is balanced by the internal force 
J-. From this we inferred the hypothesis of "weak-clustering" , namely that the typical mass 
density of dipole moments in a galaxy (after cosmological evolution) is much less than the 
baryonic density, o± <C pt>, and perhaps of the order of the mean cosmological value, <r* ~ <r*. 
Furthermore the dipolar medium is essentially static, v ~ 0. If this hypothesis is true, we 
have g ~ J- by (3.4), so the polarization field 11^ is aligned with the gravitational field g, 
i.e. the medium is polarized. Using T = £± W together with the expression of the potential 
(3.1 ), we get 

2~47rIIj_(l + 47r— j + 0(U 3 ± ) . (3.7) 



Hence the gravitational susceptibility coefficient x = A 4 ~ 1 defined by (1.2) takes the appro- 
priate form in the MOND regime, namely x(d) — ~ l + 5 , / a o + ^(5' 2 )- We conclude that (3.6) 



is equivalent to the MOND equation (1.1). (See paper I for a discussion of the Newtonian 
regime g ^> a .) Note that it is crucial that we could neglect the monopolar part of the 
dipolar medium as compared to pt>, so that galaxies appear baryonic in MOND fits of the 
rotation curves. On the other hand, the monopolar dark matter o± as we have seen plays 
the dominant role in a cosmological context. It may also help explaining the missing dark 
matter at the intermediate scale of galaxy clusters [31]. 

The weak-clustering mechanism also tells us that the evolution of the dipole moments 



should be slow. In spherical symmetry, the two last terms of (3.5) cancel each other, and 
we get d 2 U± = Ana^U^ in the MOND regime. This shows the presence of an instability, 
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with exponentially growing modes. However the unstable modes will develop on the self- 
gravitating time scale r g = a/vt/ct^, which is very long thanks to a* <C pb- Using the mean 
cosmological value a* ~ 10~ 26 kg/m 3 we get r g ~ 6 x 10 10 years. Thus this instability is not 
a problem classically. 



IV. CONCLUSION 

In conclusion, the model (i) explains the phenomenology of MOND by the physical process 
of gravitational polarization, (ii) makes a unification between the dark matter a la MOND 
and the dark energy in the form of a cosmological constant (with the interesting outcome 
that A ~ <2q), and (iii) recovers the successful standard cosmological model A-CDM at linear 
perturbation order. However the model lacks some connection to microscopic physics and 



describes the dipole moments in an effective way; notably the potential W in (3.1) is for 
the moment purely phenomenological. The model should be further tested in cosmology, 
by studying second-order cosmological perturbations where we expect a departure from A- 
CDM, by computing numerically the non-linear growth of perturbations and formation of 
large scale structures, and by investigating the intermediate scale of galaxy clusters. 



APPENDIX A: VARIATION OF THE ACTION FUNCTIONAL 

Here we provide some details on the derivation of the equations of motion and evolution 
of the dipolar fluid. They derive from an action of the general form 

ii 



S = J d 4 xv^L[J^,e,^], (Al) 

where as indicated the Lagrangian density L is a functional of the matter current J M , the 
dipole moment £ M (and its covariant derivative V^), and the covariant metric g^ u . 

We vary first the action with respect to the dipole moment Notice that in our 
Lagrangian (2.2) the dependence on V l/ ^ is only through the covariant time derivative 



In that case, denoting the conjugate momentum of the dipole by ^ = <9L/<9£ M , we obtain 
from the principle of stationary action 

BL 

^ + e ^ = ^' (A2) 

with 9 = V v u v '. Since the vector field is unconstrained, this equation is equivalent to the 
standard Lagrange equation 

„ / dL \ dL , . 

Ha^) = W (A3) 



dipolar fluid as given by (2.3). 



In the case at hands of the Lagrangian (2.2) we then obtain the equation of motion of the 



However, the variation with respect to the current is trickier because of the constraint 
that this current is conserved: V M J M = 0. We adopt a convective variational approach 
[HI HE] in which the variation SJ^ is constrained to have the form which is precisely 
induced by an infinitesimal displacement of the flow lines of J M . Denoting 5x^ the generator 
of the displacement of the flow lines we have 

5J» = bx^^ - rv u 5x» + J»V v 5x v , (A4) 
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which is automatically divergenceless: V M £ J M = 0. The variation of the action with respect 
to J M , using the fact that 8x^ is unconstrained, then yields 



^[V^-V^pJ =0 



(A5) 



with = dL/dJ^ being the momentum associated with the current. In the case of a perfect 
fluid this equation is equivalent to the usual Euler equation, where in that case is the 
current of enthalpy [47, 48J. For the Lagrangian (2.2), the equation (A5) yields the equation 



of evolution of the dipolar fluid in the form (2.4) 



Finally the stress-energy tensor is derived by variation of the Lagrangian with respect to 
the metric. We take into account the dependence of the current on the metric through 
the volume element y/—g d 4 x, which means that the so-called "coordinate" current density 
J* = \/~9J^ is the relevant metric- in dependent variable. In addition we treat the change 
in the metric that is hidden into the covariant time derivative £ M by means of the Palatini 
formula. The result is 

T »v = 2 — + gltu r L _ JP \ + u » u v ^ + v ^ L(p^)^p _ u p^m») _ £p u (^»A . (A6) 

CJQii.ij \ / \ / 



(Notice the misprint in the first dipolar term in the corresponding equation (2.22) of [36].) 
Straightforward calculations in the case of the Lagrangian density (2.2) give the explicit 



expression (2.5). 
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